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ABSTRACT 

We study population of compact X-ray sources in the Small Magellanic Cloud using 
the archival data of XMM-Newton observatory. The total area of the survey is « 1.5 
square degrees with the limiting sensitivity of ~ ICR 14 erg/s/cm 2 , corresponding to 
the luminosity of « 4.3 • 10 33 erg/s at the SMC distance. Out of ~ 150 point sources 
detected in the 2-8 keV energy band, ~ 3/4 are background CXB sources, observed 
through the SMC. Based on the properties of the optical and near-infrared coun¬ 
terparts of the detected sources we identified likely HMXB candidates, and sources, 
whose nature is uncertain, thus, providing a lower and upper limits on the luminosity 
distribution of HMXBs in the observed part of the SMC. 

The observed number of HMXBs is consistent with the prediction based on SFR 
estimates derived from the supernovae frequency and analysis of color-magnitude dia¬ 
grams of the stellar population. If, on the contrary, the true value of the SFR is better 
represented by FIR, H a and UV based estimators, then the abundance of HMXBs in 
the SMC may significantly (by a factor of as much as ~10) exceed the value derived 
for the Milky Way and other nearby galaxies. The shape of the observed distribu¬ 
tion at the bright end is consistent with the universal HMXB XLF. At the faint end, 
Lx fS 2 • 10 34 erg/sec, the upper limit on the luminosity function is consistent with 
while the lower limit is significantly flatter than the L~ oe power law. 

Key words: X-rays: galaxies - X-rays: binaries - stars: neutron - galaxies: individual: 
SMC. 


1 INTRODUCTION 

X-ray observations of the Small Magellanic Cloud (SMC) 
revealed a surprisingly rich population of High Mass X- 
ray Binaries in this close neighbour of the Mi lky Way (e.g. 
lYokogawa et al .Il2000l : lllaberl fa Pietschll2004i . At the time 
of writing, about 50 sources of this type have been iden¬ 
tified in the SMC with the majority of discovered systems 
being Be/X-ray binaries. Existence of such a rich popula¬ 
tion of HMXBs in a nearby galaxy opens unique possibility 
to study properties of popu latio n of HMX Bs. 

As has been shown bv lUrimm et ali (2003), the X-ray 
luminosity function (XLF) of HMXBs obeys, to the first 
approximation, a universal power law distribution with the 
differential slope of « 1.6, whose normalization is propor¬ 
tional to the star formation rate of the host galaxy. The 
validity of this universal HMXB XLF has been established 
in a broad range of the star formation rates and regimes 
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and in the luminosity range log(Lx) > 35.5 — 36. Based 
on the ASCA observations of the Small Magellanic Cloud 
and on the behavior of the integrated X-ray luminosity of 
distant galaxies located at redshifts z ~ 0.3 — 1. 3 observed 
b y Ch andra in the Hubble Deep Field North. fOrimin et all 
1 2003 1 tentatively suggested, that the HMXB XLF is not 
dramatically affected by the metallicity variations. 

Study of the population of high mass X-ray binaries in 
the SMC is of importance for several reasons: 

(i) The Magellanic Cloud s are known to h ave a significant 
under-abundance of metals |jWesterhindll997ft . For example, 
the interstellar medium of the SMC has a mean metallicity 
0.6 dex lower than the local ISM, the [Fe/H] of young objects 
(age< 0. 6Gyr) in the SMC is 0.5 dex lower than in the solar 
vicinity dWesterlundl Il997l l. The effects of the metallicity 
variations on the population of X-ray binaries are poorly 
understood. Study of the SMC, a galaxy with relatively well 
known chemical composition gives a unique opportunity to 
investigate such effects observationally. 

(ii) Owing to the proximity of the SMC, the weakest 
sources become reachable within a moderate observing time. 
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Figure 1. Far-infrared (IRAS, 1 00//) map of the region around 
the Small Magellanic Cloud. Circles show fields of view of the 
XMM-Newton observations used for the analysis. 


Indeed, the sensitivity of a typical Chandra or XMM- 
Newton observation, ~ 1CP 14 erg/s/cm 2 corresponds to the 
luminosity of ~ 4.3 • 10 33 erg/s at the SMC distance. This 
opens the possibility to study the low luminosity end of the 
HMXB XLF. In addition, the proximity of the SMC allows 
to study the properties of its stellar population and the star 
formation history in details. 

In the present paper we study population of X-ray bi¬ 
naries in the SMC based on the archival XMM-Newton 
data. Th e distanc e modulus of the SMC is m — M = 18.9 
llWesterlundil997l). corresponding to the distance of D ~ 60 
kpc. 

The paper is structured as follows. XMM-Newton ob¬ 
servations and data analysis are described in the section 0 
In the section [3]we discuss nature of detected X-ray sources. 
In the section Q]we describe the HMXBs search procedure 
and its results. Resulting HMXBs luminosity function is pre¬ 
sented in the section 0 The log(N)-log(S) of CXB in the 
direction of the SMC is presented in the section |fi] 


2 OBSERVATIONS AND DATA ANALYSIS 

We have selected 9 XMM-Newton archival observations with 
the pointing direction towards the SMC and with sensitivity 
better than ~ few x 10~ 14 erg/s/cm 2 in the 2-10 keV energy 
band. These observations are listed in Table □ FigQ shows 
their fields of view, overlayed on the far-infrared map (IRAS, 
100/r) of the Small Magellanic Cloud. 

The observations were processed with Science Analysis 
System (SAS) v6.0.0. After filtering out high background 
intervals we extracted images in the 2-8 keV energy band. 
This energy range was chosen to minimize the fraction of 
cataclysmic variables and foreground stars among detected 
sources, generally having softer spectra than high mass X- 
ray binaries. To improve the sensitivity of the survey, images 
from MOS1 and MOS2 detectors have been merged. If both 
MOS and PN data were available, we used the data having 


higher sensitivity. This was determined comparing the ex¬ 
posure times and taking into account that the effective area 
of PN is ss two times higher than the effective area of each 
of MOS cameras. 

2.1 Source detection 

The source detection was performed with standard SAS 
tasks eboxdetect and emldetect. The value of the threshold 
likelihood L used in emldetect task to accept or reject de¬ 
tected source was chosen as follows. We simulated a number 
of images with Poisson background (without sources). Each 
of the generated images was analyzed with the full sequence 
of the source detection procedures using different values of 
the emldetect threshold likelihood L and for each trial value 
of L the number of detected “sources” was counted. The fi¬ 
nal value of the threshold likelihood, L = 10.5, was chosen 
such that the total number of spurious detections was < 1 
per 9 images. 

The obtained images were visually inspected and the 
source lists were manually filtered of spurious sources near 
bright sources and sources whose extent undoubtedly ex¬ 
ceeded PSF size. The final merged source list contains « 150 
sources. Their flux distribution is plotted in FigED 

2.2 Energy conversion factor 

The 2-8 keV source counts were converted to the 2-10 keV 
energy flux assuming a power law spectrum with the photon 
index 1 and N/r=5T0 2 °cm~ 2 . The energy conversion fac¬ 
tors are ecfMOS=2.95-10 -11 erg/cm 2 and ecfpN=1.0-10~ n 
erg/cm 2 . We note that, as will be shown in sec. 13.11 signif¬ 
icant part of sources in our sample are background AGN. 
Therefore, in sections rm and 0 (figures 0] [U and Q 
we used energy conversion factors calculated for photon in¬ 
dex 1.7, ecfMOS=2.27-10~ n erg/cm 2 and ecfpN=8.0-10 -12 
erg/cm 2 , instead of values mentioned above. 

2.3 Boresight correction 

For several XMM observations we performed boresight cor- 
rection using optica l sources from the MCPS catalogue 
iZaritskv et, a.l ] [200211 . The correction was applied in case 
when two or more well-known sufficiently bright X-ray 
sources with optical counterparts were found. 

2.4 Correction for incompleteness 

Naturally, the point source detection sensitivity varies from 
observation to observation and, for a given observation, 
across the telescope field of view. The reasons of these vari¬ 
ations are difference in the exposure times, deterioration of 
the point spread function of the telescope with the increase 
of the off-axis angle, presence of the diffuse emission and 
background variations. These factors affect the complete¬ 
ness of the survey at low fluxes. The conventional way to 
account for these effects without setting too high conserva¬ 
tive completeness limit, is to calculate (using the sensitivity 
maps of individual observations) the survey area as a func¬ 
tion of flux. After appropriate normalisation it represents 
the fraction of sources we are able to detect as a function 
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Table 1. List of XMM-Newton observations used for analysis. 


Obs. ID 

Target 

R.A. 

J2000 

Dec. 

J2000 

Instrument 

Exposure 

ksec 

0135721501 

1ES0102-72 

16.0297 

-72.0109 

PN 

25ks 

0112880901 

CF Tuc 

13.3742 

-74.6614 

MOS1+MOS2 

80ks 

0110000201 

IKT 18 

14.9441 

-72.1583 

MOS1+MOS2 

39ks 

0110000301 

IKT 23 

16.3021 

-72.3744 

MOS1+MOS2 

58ks 

0110000101 

IKT 5 

12.3645 

-73.2215 

MOS1-I-MOS2 

53ks 

0084200101 

SMC Pointing 1 

14.1015 

-72.3590 

MOS1+MOS2 

33ks 

0084200801 

SMC Pointing 8 

13.7263 

-73.6701 

MOS1+MOS2 

42ks 

0011450101 

SMC X-l 

19.1750 

-73.4381 

MOS1+MOS2 

102ks 

0157960201 

XTE J0055-727 

13.9263 

-72.7143 

MOS2 

19ks 



Figure 2. The survey area as a function of 2-10 keV flux, calcu¬ 
lated as described in section m 


of flux, and can be used to correct the detected number of 
low-flux sources: 
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where A(S) is the survey area at flux S and Ao is the total 
(geometrical) survey area. The corrected cumulative log(N)- 
log(S) distribution can be obtained as follows: 


N(> S ) = 


E 

Sj>S 


Ao 

MSjY 


( 2 ) 


where Sj is the flux of the j-th source. 

The f lux-depende nt su r vey ar ea calculated as described 
in IShtvkovskiv fc Gilfanovl 1 2005 ) is shown in Fig. |2] From 
Fig • I^Jone can see that the incompleteness effects become im¬ 
portant at fluxes < lCP 13 erg/s/cm 2 . In the high flux limit, 
the total area of the survey equals /I t.ot — 1.48 deg 2 . The 
flux distribution, corrected for the incompleteness effects, is 
plotted as the thick histogram in FigEl 


3 NATURE OF X-RAY SOURCES IN THE 
FIELD OF THE SMC 

3.1 Background and foreground sources 

The total number of sources with the flux Fx [2 — 10 keV] > 
2.3 • 1CP 14 erg/s/cm 2 (corresponding to the luminosity 
10 34 erg/s at the SMC distance) is 151. With account for 
the survey incompleteness (Eq. 0 this number corresponds 
to N(> 2.3 • 10 -14 ) « 192 (see thick histogram in Fig. 0. 
According to the CXB log(A) — log(S') determined by 
iMoretti et al (2003 ). the total number of CXB sources ex¬ 
pected in the held of 1.48 deg 2 is Ncxb ~ 148. From the 
comparison of these numbers it is obvious that about 3/4 of 
sources in our sample are background AGNs. 

A significantly less important source of contamination is 
X-ray sources associated with foreground stars in the Galaxy 
(no known Galactic X-ray binaries were located in the held 
of view of XMM observations). Foreground stars are re¬ 
moved by our filtering procedure, as described in section 0 


3.2 Low mass X-ray binaries 

Given the limiting sensitivity of the survey, Fx ~ (1 — 3) ■ 
10“ 14 erg/s/cm 2 , corresponding to the luminosity Lx ~ 4.3- 
10 33 — 1.3 ■ 10 34 erg/s at the SMC distance, the intrinsic 
SMC sources are dominated by X-ray binaries. Their total 
number is proportional to the stellar mass (LMXB) and star 
formation rate (HMXB) of the galaxy. 

The total stellar mass of the SMC can be estimated 
from the i ntegrate d optica l lumi nosi ty. According to RC3 
catalog ide Vaucouleurs et alJ I 199 1). the reddening cor¬ 
rected V-band magnitude of the SMC equals Vto ~ 1.92, 
corresponding to the total V-band luminosity of Ly « 
4.6 • 10 s Lc tv From t he dereddened optical color (B — 
V )to ~ 0 -36 ide Vaucoule urs et alJI 199 if) and using results 
of lBell fc de Jon^ ~ l20^1) ~the V-band mass-to-light ratio is 
(M/L)v ~ 0.59 in solar units, giving the total stellar mass 
of the SMC M* ~ 2.7 ■ 10 8 Mq. For this stellar mass, using 
results of iGilfanovI 1 2004 ). about one LMXB with luminos¬ 
ity Lx > 10 35 erg/s is expected in the entire galaxy. Such a 
prediction agrees with observations as no LMXB candidates 
are known in the SMC at present. 
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Figure 3. Cumulative log(N)-log(S) distribution of detected 
point-like sources, excluding several known foreground stars and 
rotation powered pulsars (section [3J. The thin and thick his¬ 
tograms show, respectively, the observed distribution and the dis¬ 
tribution corrected for the incompleteness effects, as described in 
section 12.41 The thick grey line sho ws the log(N)-log(S) of CXB 
sources according to iMore tti et aj £ 2003 1. The brightest source 
SMC X-l with flux Sx = 6.15 • 10 -10 erg/s/cm 2 is situated far 
outside the flux range of the plot. 

3.3 Star formation rate in the SMC 

Determination of a star formation rate in the SMC, a diffuse 
galaxy with low opacity is a nontrivial task. 

Below we compare estimates of the star formation rate 
in the SMC, obtained from different SFR indicators. The 
SFR. values refer to the star formation rate of stars in the 
0.1-100 Mg range, assuming Salpeter’s IMF. The discussion 
of results of diff erent star formation indicators could be also 
found in IWilke et_aljjj 200^W see also di s cussio n of the SFR 
in t he LMC bvIShtvkovsk R fc Gilfanov! ll2005ll 1 . 

IWilke et all llgOO^Cbased on ISOPHOT observations 
and calibration of lKennicuttl (IlDDSTl obtained th e star forma- 
tion rate in the SMC SFR(IR)=0.015 /vr. [Wilke et alJ 
i 2004 1. however, noted that this value underestimates the 
SFR of the SMC because some stellar radiation, due to low 
opacity of the SMC, can escape the star formation regions 
directly, and gave a more realistic estimate of 

SFR(IR) = 0.05 Mq/ yr. (3) 

With the total infrared luminosity of the SMC lWilke et ahl 
J‘2004T L : n = 3.21-10 41 erg/s and using the equation (5) from 
iBelll f: 200.'fl which accounts for the uncertainty mentioned 
above, we obtained 

SFR(IR) = 0.044 Mq/ yr. (4) 

Note, that this calibration should be applie d extr emely cau¬ 
tious for such a low-luminosity galaxy feell200 3.!. 

The H a luminosity-based indicators can also underesti¬ 


mate the star formation rate in the SMC due to the reasons 
mentioned above. Based on H a luminosity of the SMC and 
applying extinction correction FKennicutt et all ( 1995 1 esti¬ 
mated the star formation rate of 

SFR(H q ) = 0.044 M e / yr. (5) 

iKennicuttl i lODll i gave an conservative upper limit of 

SFR(H a ) < 0.1 Mq/ yr. (6) 

An SFR estimator relatively free of this uncertainty 
is the one based on combined UV and FIR emission. 
With the UV flux measu red by the D 2B-Aura satellite 
liVangioni-Flam et alll98flil at A = 1690A corrected for fore¬ 
ground extincti on and unc orrect ed ISOPHOT FIR flux we 
obtain with the lKennicutt 01998 4 calibration 

SFR(UV + FIR) = 0.063 Mq/ yr. (7) 

This value is consistent with the SFR(UV) derived from 
fl ux corrected for bo th for eground and internal extinction 
iVangioni-Flam et al.lll980 ): 

SFR(UV, ext.corr.) = 0.063 Mg/j/r (8) 

iFilipovic et alJ dl99Si) . from comparison of discrete ra¬ 
dio (Parkes telescope) and X-ray (ROSAT) sources, esti¬ 
mated a number of SNRs in the SMC (14). From the age 
- radio flux density relation they estimated an SNR birth 
rate of one SNR in 350 ± 70 yrs. We convert this value to 
the SFR in the mass range 0.1-100Mq, assuming Salpeter’s 
IMF and taking into account that only stars with M^ 8 Mq 
can result in a radio SNR: 

SFR(SNR birth rate) « 0.38±g;^ Mq/ yr. (9) 

_ Similar value has been obtained by lHarris fc Zaritskvl 

(2004) who studied star formation history of the SMC fit¬ 
ting the color-magnitude diagrams of stellar population with 
theoretical isochrones: 

SFR(CMD) « 0.37 Mq/ yr. (10) 

This number corresponds to recent star formation which 
took place ~ 4 — 5Myr ago. 

Thus, the SFR values derived using different star for¬ 
mation rate proxies can differ by the factor of as much as 
~ten. The “standart” indicators (IR, UV, H a ) result in the 
rather low value of ~ 0.05 Mg/yr, whereas the estimates 
based on the radio supernova rate and the analysis of the 
color-magnitude diagram point towards much larger star for¬ 
mation rates, in the ~ 0.3 —0.4 Mg/yr range. We note, that 
the methods based on the analysis of the color-magnitude 
diagrams are generally believed to be the most reliable and 
accurate. As the controversy between different SFR indica¬ 
tors can not be resolved at present, we shall assume in the 
following that the true value of SFR is somewhere in the 
range: 

SFR(SMC) « 0.05 - 0.4 Mq/ yr. (11) 

3.4 High mass X-ray binaries 

For the star formation rate of SFR=0.05 — 0.4 Mq yr -1 , 
using HMXB-SFR relation of lGrimm et, alJ (l2003ll (see also 
the discussion of th e normalizat ion in the HMXB-SFR cal¬ 
ibration by IShtvkovskiv fc Gilfanovll2005 1 we predict from 
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~ 6 to ~ 49 HMXBs with luminosities > 10 35 erg s _1 in 
the whole SMC. This confirms that the population of X-ray 
binaries in the SMC is dominated by HMXBs. 

To predict the number of HMXBs in our sample of X- 
ray sources, we estimate the star formation rate in the part 
of the SMC covered by XMM pointings from IRAS in f rared 
maps provided by Sky View llMcClvnn. Scollick fc Whitel 
119961 1. Calculating the far infrared flux according to the 
formula FIR= 1.26 ■ 10 -11 (2.58S , 6o ( j + Sioof), where FIR 
is flu x in erg/s/cm 2 and S\ is flux in Jy llHelou et alJ 
Il985tl . and integrating the IRAS maps we obtained the total 
SFR= 0.011MQ/yr. This number is a factor of ~ 4.5 - 7 - 36.4 
smaller than eq. im because of the uncertainties mentioned 
in sec. roi To make it consistent with our determination of 
the SFR in the SMC (eq 1111 . we simply multiply the FIR- 
based values by the correction factor of fa 4.5 -r 36.4. The 
thus calculated star formation rate of the part of the SMC 
covered by XMM observations, is: 

SFR(XMM) fa 0.018 4- O.15M 0 /yr. (12) 

With the above value of SFR we predict about ~ 2 — 18 
HMXBs with luminosities ^ 10 35 erg s -1 in the observed 
part of the SMC. The error in this number accounts for 
uncertainty of the SFR estimate and does not include the 
Poisson fluctuations. 

3.5 The Log(N)—Log(S) distribution 

As demonstrated above, the population of compact X-ray 
sources in the SMC field is dominated by two types of 
sources - background AGN and high mass X-ray binaries in 
the SMC. Their log(N) log(S) distributions in the flux range 
of interest can be described by a power law with the differen¬ 
tial slopes fa 2.5 (CXB) and fa 1.6 (HMXBs). Due to signifi¬ 
cant difference in the slopes, their relative contributions de¬ 
pend strongly on the flux. At large fluxes, Fx > (2—3)-10“ 13 
erg/s/cm 2 [Lx > 10 35 erg/s) the X-ray binaries in the SMC 
prevail. On the contrary, in the low flux limit, e.g. near the 
sensitivity limit of our survey, Fx ~ 10~ 14 erg/s/cm 2 , the 
majority (~ 3/4) of the X-ray sources are background AGN. 

This is illustrated by Fig. [3 showing the observed and 
corrected for incompleteness log(N)-log(S) distribution of 
all sources from the final source list. The corrected log(N)- 
log(S) distribution agrees at low fluxes with that of CXB 
sources. At high fluxes, there is an apparent excess of sources 
above the numbers predicted by the CXB sources log(N)- 
log(S), due to the contribution of HMXBs. 


4 IDENTIFICATION OF HMXB CANDIDATES 

To filter out contaminating background and foreground 
sources, we use t h e same filtering procedure as described by 
IShtvkovskiv V Gilfa.uovl (' 2005 1. It is based on the fact that 
optical emission from HMXBs is dominated by their optical 
companions, whose properties, such as absolute magnitudes 
and intrinsic colors are sufficiently well known. 

4.1 Optical properties of HMXBs counterparts 

High mass X-ray binaries are powered by accretion of mass 
lost from the massive early-type optical companion. The 


mechanism of accretion could be connected either to (i) 
strong stellar wind from an OB supergiant (or bright giant) 
or (ii) eq uato rial circum stellar di sk arou nd Oe or Be typ e 
star (e.g. ICorbetl Il986l : Ivan Paradiis fc McClintockl Il994l) . 
In the case of Small Magellanic Cloud where the only su¬ 
pergiant system SMC X-l is known at present time, the 
HMXB s with Be c ompanion seem to give the main contri¬ 
bution ( Coe et a l. i200 4l. 

Given the large number of well known high mass X-ray 
binaries in the SMC, the selection criteria for the optical 
counterparts of such sy stems can be determined straight¬ 
forwardly. For example, ICoe et ah (12004 4 studied the opti¬ 
cal properties of Be/X-ray pulsars in the SMC. The major¬ 
ity of HMXBs with identihed optical companions in their 
sample (32 out of 34) have apparent V-band magnitudes 
my = 14 — 17. The exceptions are supergiant system SMC 
X-l with my = 13.2 and X-ray pulsar CXOU J010042.8- 
721132 with my = 18.01 (OGLE catalogue [Udffisk^^ffi. 

which has been previously classified hv TLmn^lcah 
(1 20021) as a possible anomalous X-ray pulsar (AXP). 

This picture is consistent with the expectations based 
on the positions of possible optical counterparts on the 
Hertzsprung-Russel dia gram , distance modulus of the SMC, 
(m — M) 0 = 18.9 dWesterlundll997i) . and foreground and in¬ 
trinsic ex tinction toward s the SMC, Ay ~ 0.46 (for hot pop¬ 
ulation) dZaritskv et ali 1*200 2). OB supergiants and bright 
giants (luminosity classes I—II) have absolute magnitudes 

My ~ — 7-4, corresponding to apparent magnitudes in 

the range my ~ 12.5 — 15.5. The position of Be stars in 
the H R diagram is close to the main sequence. Apparent 
magnitudes my < 17 — 18 correspond to spectral classes 
earlier than B3-B5 for main sequence stars (B5-B7 for gi¬ 
ants). Thus we can conclude, that HMXBs in the SMC with 
main sequence Be optical companions have spectral classes 
earlier than B3-B5 (B5-B7 for giants), in agreement with 
expectations based on the optical properties of HMXBs in 
our Galaxy. 

As potential optical counterparts of HMXBs belong to 
spectral classes earlier than ~ B3-B5, their intrinsic optical 
and near-infrared colors are constrained by B — V < —0.20, 
J — K < —0.16. With account for the interstellar redden¬ 
ing, both apparent colors are expected to be < 0T_ — 0.2. 
This i s consistent with B-V colors of HMXBs from lCoe et alJ 
t:2004h sample. In infrared part of the spectrum, however, 
Be stars are known to exhibit an infrared continuum excess, 
caused by free-free and free-bound emission within the disk 
dCehrz et alll974lh Therefore, J-K colors of Be stars can be 
higher than values qouted above. To avoid this uncertainty, 
in identification procedure we use mainly B-V colors, ap¬ 
plying J-K identification criteria only in extreme cases (in 
fact we used this criteria alone only for one source). 

4.2 Catalogs and selection criteria 

We used the following optical and near-infrared catalogs: 

(i) Ma gellanic Clouds Photometric Survey: the SMC 
(MCPS) dZaritskv et alJl200 2) 

(ii) Gu ide Star Catalog, version 2.2.1 (GSC2.2.1) 

dMorrison fc McLeanll200ll f 

(iii ) The C CD survey of the Magellanic Clouds 
llMassev et aljEoO^l 
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(iv) 2-micron All Sky Survey (2MASS) llCutri et all2Qo3l 

(v) Emission-line sta rs and PNe in the SMC 
jMevssonnier et al.lftogi^l 

Cvi) A 2d F su rvey of the Small Magellanic Cloud 
llEvans et alJl2004h 

As a first step, we cross-correlated the XM M-Newton X- 
ray source list with MCPS, CSC and lMassev et all ll2002tl 
optical catalogs, using search radius 4" and the following se¬ 
lection criteria: 12.0 < V < 18.0, B — V < 0.6. If optic a l ob¬ 
ject w as present in MCPS and either CSC or lMassev et alJ 
(' 2002 1 catalogs, the preference was given to MCPS, as it 
has a higher photometric accuracy (especially at low mag¬ 
nitudes). The color limits are higher than possible colors of 
HMXBs optical companions (B-V<0.1-0.2) and were chosen 
to account for limited photometric accuracy of optical cata¬ 
logues. We found optical counterparts for 54 X-ray sources 
using these selection criteria. 

On the second stage, we cross-correlated all X-ray 
sources having optical counterparts with near-in frared cata- 
lo g 2MASS, ca t alog o f lMevssonnier et al . 1 199311 and catalog 
of lEvans et alJ ( 200411 using the same search radius as before. 

Then the following filtering procedure was applied. 

(i) In case when source was present only in the GSC cat¬ 
alogue (which has rather low photometric accuracy), we ap¬ 
plied filtering criteria based on the near infrared data. In 
fact, there was only one source with J-K«1.2 rejected by 
this criteria. 

(ii) All X-ray sources with low X-ray-to-optical flux ra¬ 
tio Fx/F 0 pt < 10 -3 were rejected. The optical flux was 
calculated from the V-band magnitude, F op t = 8.0 • 10 -6 ■ 
10 -mv / 2 ' 5 erg/s/cm 2 . Such low Fx/F 0pt < 10 -3 ratios are 
typical for foreground stars but not for X-ray binaries. All 
confirmed HMXBs in our sample have Fx/F op t > 3 ■ 1CP 3 
among which > 90% have Fx/F op t > 10~ 2 . 

4.3 Search radius 

The choice of the search radius plays a crucial role in the 
procedure of search for potential HMXB candidates based 
on their optical companions. Too small value of the search 
radius may lead to incompleteness of resulting HMXBs list 
(especially at low luminosities), while too large value may 
lead to a significant fraction of chance coincidences. Given 
the spatial density of optical stars (whose optical colors and 
magnitudes satisfy our selection criteria) in the MCPS cata¬ 
logue, the number of detected X-ray sources and the search 
radius of 4”, we estimate the number of chance coincidences 
to be « 14. On the other hand, a typical error on the po¬ 
sition of a low-luminosity source is a ~ 2" (lcr confidence). 
Therefore, in ideal situation with the search radius of 4" we 
miss only ~ 13% of low-luminosity sources. The real fraction 
of missed optical counterparts, however, could be somewhat 
higher as it was obtained assuming ideal astrometric accu¬ 
racy. Our choice of the search radius is further justified in 
the Appendix. 

4.4 Identification results 

After the filtering, 50 X-ray sources remained in the list 
of potential HMXB candidates (Table [5J. As follows from 


sec 14. 3 I « 14 of these sources can be chance coincidences. 
Therefore, we further split the filtered list into reliable 
HMXBs and the sources whose nature is less certain. We 
include in the list of reliable HMXB candidates only sources 
which have either detected X-ray pulsations or optical emis¬ 
sion lines, or sources whose spectral class is well-determined 
based on their optical spectrum. The latter two ide ntifica¬ 
tion criteria are mainly based on iMevssonnier et alJ (Il993h 
and lEvans et, al.l (1200411 catalogues. 

There were 32 sources satisfying at least one of these 
criteria. Remaining 18 sources are considered to be sources 
of uncertain nature. 

The resulting likely HMXBs candidates and sources of 
uncertain nature are listed in Table [5] All known HMXBs 
located in the field of view of the XMM-Newton observations 
are among these sources. 

Comments on the individual sources: 

#1: To minimize the influence of pile-up on the flux 
of SMC X-l, for this source we used data from EPIC PN 
camera operating in Small Window mode. 

#17 = RXJ0049.2-7311: Optical colors consistent with 
HMXB n ature. This source was identified with an bright H Q 
object bv lGoe et al.l t.200j) and is probably associated with 
nearby X-ray pulsar AXJ0049-732 (p=9s). 

#34: Optical colors consistent with HMXB nature, how¬ 
ever, infrared colors are relatively high, J-K=0.97. 

#35: O ptic al colors consistent with HMXB nature. 
lLamb et alJ (1200 24 classified this source as a possible anoma¬ 
lous X-ray pulsar (p=8s). 

#38: This source is situated inside the NGC 330 stellar 
cluster. As a stellar density is high there, we consider the 
photometry to be unreliable. 

#49: Optical colors consistent with HMXB nature. 
This source is situated in 24" from source AXJ0048.2-7309. 
lYokogawa et al • l2003h based on the presence of emission- 
line object in lMevssormier et, alJ lll993ll catalogue within the 
ASCA error circle of 40" (90% confidence) and X-ray spec¬ 
trum classified this ASCA source as a possible Be X-ray 
pulsar. The posititon of th e XMM source is i nconsi stent 
with emission-line object in IMevssonnier et al . Il99 3'f and 
therefore either optical identification of ASCA source was 
incorrect either XMM source is unrelated to ASCA source. 
Therefore we consider this source as a source of uncertain 
nature. 


4.5 Completeness 

The completeness of the list of HMXB candidates is defined 
by the following factors: 

(i) The completeness of the optical catalogues. 
The initial search for optical counterparts is based 
on MCPS and GSC2.2 catalogs. In general, the 
MCPS catalog u e is thought to be complete for V<20 
iZaritskv et al, l2002ll . The GSC2.2 is a magnitude 
selected (V # 19.5) subset of the GSC-II catalog 

(http://www-gsss.stsci.edu/gsc/gsc2/GSC2home.htm). 

The lat ter is com plete to J = 21 at high galactic latitudes 
jMorrison fc McLeanfeoQll l. 

Completeness of both catalogs is known to break down in 
the crowded regions. As no sensitivity maps for the optical 
catalogs exist, a quantitative estimate of the completeness 


































High mass X-ray binaries in the SMC: the luminosity function 7 


Table 2. List of HMXB candidates 


# 


R.A. 

Dec. 

Lx 

(1) 

my 

B-V 

Fir/F opt 

id 

Comments 








erg 

s" 1 






Likely HMXB candidates 

i 

01 

17 

05.2 

-73 

26 

38 

3.8 • 

10 38 

13.2 

-0.14 

0.55 

p, ob 

SMC X-l P=0.717s (2) 

2 

00 

47 

23.3 

-73 

12 

28 

1.4 ■ 

10 36 

16.0 

0.08 

0.80 

p 

RX(AX)J0047.3-7312 P=263s (2) 

3 

00 

51 

52.3 

-73 

10 

34 

1.1 ■ 

10 36 

14.5 

-0.07 

0.14 

p,em 

RXJ0051.9-7311=AXJ0051.6-7311 P=172s (2) 

4 

00 

49 

42.0 

-73 

23 

14 

6.4 ■ 

10 35 

14.9 

0.15 

0.12 

p,em 

RXJ0049.7-7323=AXJ0049.5-7323 P=755.5s (2) 

5 

00 

57 

50.4 

-72 

07 

55 

6.1 ■ 

10 35 

15.7 

-0.05 

0.26 

p,em 

CXOUJ005750.3-720756 P=152.3s (2) 

6 

00 

54 

55.9 

-72 

45 

10 

5.2 ■ 

10 35 

15.0 

-0.03 

0.14 

p,em 

AXJ0054.8-7244 P=500s(2,3) 

7 

00 

57 

49.6 

-72 

02 

35 

4.6 ■ 

10 35 

15.7 

-0.12 

0.19 

p,em 

RXJ0057.8-7202=AXJ0058-72.0 P=280.4s (2) 

8 

01 

01 

20.8 

-72 

11 

17 

4.0 ■ 

10 35 

15.6 

-0.19 

0.15 

p,em 

RXJ0101.3-7211 P=455s (2) 

9 

01 

01 

52.1 

-72 

23 

34 

4.0 ■ 

10 35 

14.9 

-0.01 

0.084 

em, ob 

AX J0101.8-7223(2) 

10 

00 

55 

18.5 

-72 

38 

52 

3.6 ■ 

10 35 

15.9 

0.15 

0.18 

p 

P=702s (3) 

11 

00 

59 

20.9 

-72 

23 

15 

3.5 ■ 

10 35 

15.0 

-0.01 

0.075 

p 

XMMUJ005921.0-722317 P=202s (4) 

12 

00 

50 

57.3 

-73 

10 

09 

3.3 • 

10 35 

14.4 

0.08 

0.040 

em 

RXJ0050.9-7310=AXJ0050.8-7310 (2) 

13 

00 

50 

44.6 

-73 

16 

05 

2.9 ■ 

10 35 

15.5 

-0.11 

0.097 

p,em 

RXJ0050.8-7316 AXJ0051-733 P=323.2 (2) 

14 

01 

03 

13.4 

-72 

09 

13 

2.7 ■ 

10 35 

14.8 

-0.06 

0.053 

p,em 

SAXJ0103.2-7209=AXJ0103.2-7209 P=345.2 (2) 

15 

00 

58 

12.4 

-72 

30 

50 

2.0 ■ 

10 35 

14.9 

0.11 

0.038 

em 5 

RXJ0058.2-7231 (2,5) 

16 

00 

55 

28.5 

-72 

10 

57 

1.8 ■ 

10 35 

16.8 

-0.12 

0.21 

p 

CXOUJ005527.9-721058 P=34s (6) 

17 

00 

49 

13.7 

-73 

11 

38 

1.3 ■ 

10 35 

16.4 

0.19 

0.10 

em 

= (?)AXJ0049-732(P=9s) RXJ0049.2-7311 (2,7,10) 

18 

00 

48 

34.3 

-73 

02 

31 

9.9 ■ 

10 34 

14.8 

0.0 

0.023 

em, ob 

RXJ0048.5-7302 (2) 

19 

00 

54 

56.4 

-72 

26 

48 

8.1 ■ 

10 34 

15.3 

-0.05 

0.023 

p,em 

XTEJ0055-724 P=59s (2) 

20 

01 

05 

55.3 

-72 

03 

51 

6.6 ■ 

10 34 

15.7 

-0.06 

0.028 

em 

RXJ0105.9-7203 AXJ0105.8-7203 (2) 

21 

01 

01 

37.4 

-72 

04 

18 

6.1 ■ 

10 34 

16.3 

-0.16 

0.044 

em 

RXJ0101.6-7204 (2) 

22 

00 

56 

15.2 

-72 

37 

54 

4.9 ■ 

10 34 

14.6 

0.13 

0.0075 

em, ob 


23 

00 

57 

36.1 

-72 

19 

33 

4.8 ■ 

10 34 

16.0 

0.01 

0.027 

em,p 

CXOUJ005736.2-721934 P=565s (2) 

24 

00 

57 

24.0 

-72 

23 

57 

4.0 ■ 

10 34 

14.7 

-0.07 

0.0068 

ob 


25 

00 

49 

29.9 

-73 

10 

58 

3.5 ■ 

10 34 

16.2 

0.2 

0.023 

em 

= (?)AXJ0049-732(P=9s)(2,7) 

26 

01 

05 

08.1 

-72 

11 

48 

3.5 ■ 

10 34 

15.7 

-0.08 

0.016 

em 

RXJ0105.1-7211 AXJ0105-722 P=3.34s? (2) 

27 

01 

03 

37.7 

-72 

01 

35 

3.4 ■ 

10 34 

14.7 

-0.11 

0.0055 

em 

RXJ0103.6-7201 (2) 

28 

01 

19 

39.0 

-73 

30 

14 

3.3 • 

10 34 

15.9 

-0.07 

0.016 

em 

(4) 

29 

01 

01 

03.0 

-72 

07 

01 

2.7 ■ 

10 34 

15.8 

-0.07 

0.013 

em 

RXJ0101.0-7206 P=304.5 (2) 

30 

01 

00 

30.3 

-72 

20 

32 

2.1 ■ 

10 34 

14.6 

-0.06 

0.0034 

em, ob 

XMMUJ010030.2-722035 (2) 

31 

00 

56 

05.8 

-72 

21 

57 

2.0 ■ 

10 34 

15.9 

-0.04 

0.097 

p,em 

XMMUJ005605.2-722200 P=140.1 (2) 

32 

00 

50 

47.9 

-73 

18 

12 

8.7 ■ 

10 33 

15.6 

0.12 

0.0035 

em 


Sources of uncertain nature 

33 

00 

55 

34.9 

-72 

29 

05 

2.1 ■ 

10 35 

14.7 

-0.04 

0.031 



34 

00 

48 

18.8 

-73 

21 

00 

1.1 ■ 

10 35 

16.2 

0.25 

0.069 



35 

01 

00 

42.9 

-72 

11 

32 

1.1 ■ 

10 35 

17.7 

-0.03 

0.28 

p 

CXOUJ010043.1-721134 AXP? P=8s (8) 

36 

00 

54 

32.2 

-72 

18 

09 

7.2 ■ 

10 34 

16.6 

-0.08 

0.067 



37 

00 

54 

03.6 

-72 

26 

30 

6.6 ■ 

10 34 

14.9 

0.01 

0.014 



38 

00 

56 

18.7 

-72 

28 

02 

3.9 ■ 

10 34 

15.3 

-0.40 

0.012 


NGC330 stellar cluster 

39 

01 

00 

37.2 

-72 

13 

16 

2.5 ■ 

10 34 

16.7 

-0.17 

0.026 



40 

00 

48 

33.4 

-73 

23 

55 

2.5 ■ 

10 34 

17.7 

0.07 

0.064 



41 

00 

54 

41.1 

-72 

17 

20 

2.1 ■ 

10 34 

15.7 

-0.11 

0.0090 



42 

00 

47 

19.9 

-73 

08 

22 

2.0 ■ 

10 34 

17.3 

0.08 

0.036 



43 

00 

53 

35.8 

-72 

34 

24 

1.7 ■ 

10 34 

17.2 

-0.07 

0.029 



44 

01 

04 

37.4 

-72 

06 

30 

1.4 ■ 

10 34 

18.0 

0.13 

0.048 



45 

01 

00 

22.9 

-72 

11 

27 

1.2 ■ 

10 34 

17.0 

-0.13 

0.018 



46 

00 

48 

04.0 

-73 

17 

01 

1.2 ■ 

10 34 

17.0 

0.21 

0.016 



47 

00 

49 

05.4 

-73 

14 

10 

1.0 ■ 

10 34 

18.0 

0.03 

0.034 



48 

01 

03 

28.3 

-72 

06 

51 

1.0 ■ 

10 34 

16.5 

-0.15 

0.0086 



49 

00 

48 

14.9 

-73 

10 

03 

8.6 ■ 

10 33 

15.3 

0.26 

0.0025 


=(?)AXJ0048.2-7309 (2,9) 

50 

00 

50 

20.2 

-73 

11 

17 

8.6 ■ 

10 33 

17.0 

0.11 

0.012 




(1) - 2—10 keV band, a ssumin g distance of 60 kp c; (2) — seejHaberl^^Pietsch| | |20C)4|) and references in side; (3)_jHaber^^aL|Oo04); 
mjMindir^]j20p; (5) -fedgej^Copj { 20031 ; (6) - lEdge 6*^^11720041) : (7) — iFilipovic et al . ! 2003 ); (8) —|Lam^^^l!fj200^7(9) — 
lYokogaw^^t^lTl2Q0^ : (1 0) -|Co^ Co mments: p - source shows pulsations; em - source has Baimer emission lines in its 
optical spectrum (based on I M e vssonnie^e^anTlO 9 3h catalogue, unless otherwise is noted); ob — source is present in lEvans et al 1 <2004ft 

catalogue and belongs to OB spectral class; 
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Figure 4. The incompleteness-corrected XLFs of HMXB can¬ 
didates in the SMC obtained with two different HMXB search 
methods (see sec[4hj. Two solid histograms correspond to upper 
and lower limits on XLF obtained with search method based on 
the optical properties of HMXBs. Two dashed histograms corre¬ 
spond to upper and lower limits on XLF obtained with search 
method based on the X-ray hardness ratio. It is clear from this 
figure, that both methods give similar result for sources with lu¬ 
minosities Lx > 2 — 3 • 10 34 erg/s. The brightest source SMC 
X-l with luminosity L= 3.8 ■ 10 38 erg/s is situated outside the 
luminosity range of the graph. 


of the initial counterpart search is impossible. However, the 
quoted completeness limits of both catalogs are 2 — 3 mag 
better than chosen threshold of 18 mag for the optical coun¬ 
terparts search. This suggests that the completeness of the 
optical catalogs is unlikely to be the primary limiting factor. 

The MCPS catalogue also has “stripe-like” gaps with 
zero stellar density. We checked that the majority of X-ray 
sources do not coincide with this “gaps”. For few sources 
which were close or inside the gaps w e used two additional 
catalogs: OGLE lllldalski et, alj [lji9j^_and (ifCXfLE had no 
data in this region! IJSNO-B~ ilMonet et all2003lb 

(ii) The efficiency of the initial search due to statistical 
and systematic uncertainties in the positions of X-ray and 
optical sources. This is probably one of the major limiting 
factors. As has been noted in sec roi the search radius of 
4" assuming ideal astrometric quality corresponds to detec¬ 
tion of ~ 87% of the weakest sources, while for moderate 
luminosity and bright sources this fraction would be close 
to 100%. In Appendix we show, that in real case (i.e. with 
all uncertainties taken into account), with the chosen value 
of the search radius we miss not more than few sources. 

(iii) The filtering procedure applied to the optical 
matches found in the initial search. This procedure is based 
mainly on observed optical and near-infrared properties of 
known HMXBs in the SMC and Milky Way. The distri¬ 
bution of V-band magnitudes and B-V optical colors of 


Figure 5. The upper and lower limits on the XLF of HMXBs 
in the SMC. The solid histogram shows XMM-Newton results 
scaled up by the factor of SFRtot/SFRxMM ~ 2.7 according 
to the contribution of the XMM-Newton survey area to the to¬ 
tal SFR of the SMC. The dashed histogram is the ASCA XLF 
from lYokogawa et all l2003tb The shaded area shows the lumi- 
n osity distribution pr edicted from the “universal” HMXB XLF 
of lGrimm et al.l J2QQ3I) extrapolated towards low luminosities. Its 
width reflects the uncertainty of our knowledge of the SFR in the 
SMC, as discussed in the section f3.3l 

HMXBs in the SMC have maxima far inside our filtering lim¬ 
its (see Fig. [HJ, meaning that we detect the vast majority of 
HMXBs unless some significantly different (and previously 
undetected) type of these systems exists. 

Thus, we estimate the completeness of our HMXBs sam¬ 
ple to be very high with the uncertainties in the positions of 
X-ray sources being the main factor of incompleteness. We 
note, that as this factor is flux-dependent, it affects most 
strongly the low-luminosity part of our sample. 

4.6 Comparison with other HMXB identification 
methods 

The identification method based on the properties of opti¬ 
cal counterparts is not the only way to identify the popula¬ 
tion of HMX Bs in the background of c ontaminating sources. 
For example, lYokogawa et al.| ( 2000 ! found that hardness 
ratios (HRs) of all known X-ray binary pulsars (XBPs) 
in the Magellanic Clouds fall into narrow “XBP region”, 
HRxbp(ASCA)= 0.2 — 0.6. In order to compare these two 
methods, we apply selection criteria based on hardness ratio 
to our sample of X-ray sources. 

We use XBPs selection cri teria similar to those de¬ 
scribed in lYokogawa et ah] 12000). The hardness ratio is de¬ 
fined as HR = (H — S)/(H + S), where H and S are count 
rates in the 2-7 keV and 0.7-2 keV energy bands respec¬ 
tively. All X-ray sources with HRs corresponding to photon 
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ind exes a = 0.4 — 1.5 (corresponding “XBP region” defined 
by lYokogawa et alj ( 200Clh ) are classified as X-ray binary 
pulsar candidates. X-ray sources which have slightly softer 
spectra, a = 1.5 —2.0, are considered as sources of uncertain 
nature. Higher photon indexes are untypical for HMXBs, es¬ 
pecially for Be/X-ray binaries which seem to constitute the 
majority of HMXBs population in the SMC. Luminosity dis¬ 
tributions obtained using these two methods are shown in 

Fig. El 

As is clear from Fig.E] the results of both methods agree 
for sources with Lx > 2 — 3 • 10 35 erg/s, while for the low- 
luminosity sources the method based on hardness ratio gives 
significantly higher number of sources. This could be a re¬ 
sult of large statistical uncertainty in the hardness ratios for 
low flux sources, leading to a “leakage” of background AGN 
to the XBP region. In addition, the population of highly ab¬ 
sorbed AGN with truly hard spectra becomes important at 
low fluxes and compromises the simple hardness ratio based 
selection criteria. 

5 POPULATION OF HMXBS IN THE SMC 

5.1 The luminosity function of HMXBs in the 
SMC 

The incompleteness-corrected luminosity distribution of 
HMXBs in the SMC is shown in Fig. |S] In order to repre¬ 
sent the HMXB population of the entire galaxy, the observed 
XLF has been scaled up according to the relative contribu¬ 
tion of the area covered by the XMM-Newton survey to the 
total SFR of the SMC. The upper and lower histograms be¬ 
low ~ 2 • 10 35 erg/sec correspond to all sources from TablcEl 
and to the likely HMXB candidates respectively. These two 
histograms provide upper and lower limits on the true X-ray 
luminosity function of HMXBs in the observed part of the 
SMC. 

Fig. El also show s the XLF of HMXB candidates in 
the SMC obtained bv lYokogawa et’ahl feoQft) based on the 
ASCA data. The ASCA XLF refers to the entire SMC. 
We note that removal of the sources located outside XMM 
pointings does not change its shape significantly. The up¬ 
per ASCA histogram in the figure corresponds to all sources 
from lYokogawa et alJ ll2003ll excluding SNR, AGN and fore¬ 
ground stars - i.e. represents all sources which proper¬ 
ties do not contradict to HMXB nature. The lower his¬ 
togram shows the luminosity distribution of likely HMXB 
candidates - confirmed X-ray binary pulsars and candidates 
and confirmed n o n-pu lsating HMXBs and candidates (see 
lYokogawa et al ] 2003j for details). As is evident from the 
figure, the XMM-Newton and ASCA XLFs are consistent 
with each other. The K-S test for two distribution gives the 
probability of ~ 50 — 70%. 

5.2 Abundance of HMXBs in the SMC 

The s haded area in Fig.EJshows the universal HMXB XLF of 
I Grimm et al . (2003) and its extrapolation towards low lumi¬ 
nosities. The width of the shaded area reflects uncertainties 
in our knowledge of the star formation rate in the SMC as 
discussed in the section roi As is clear from Fig. El the ob¬ 
served abundance of HMXBs in the SMC is consistent with 



Figure 6. Distribution of the V-band magnitudes of the HMXB 
candidates in the SMC. The thick histogram corresponds to the 
likely HMXB candidates from the upper part of the table E] The 
thin histogram corresponds to the uncertain sources from the 
lower part of the table El 

the upper range of the predicted numbers, which is based 
on the SFR estimates derived from the supernovae rate and 
analysis of the stellar color-magnitude diagram. If, on the 
contrary, the true value of the SFR is better represented by 
FIR, H a and UV based estimators, then the abundance of 
HMXBs in the SMC may exceed significantly, by a factor 
as much as ~ 10, the va lue derived for the M ilky Way and 
other nearby galaxies bv lGrimm et al ] (i2003i) . Note_that_J,he 
analys is and interpretation of the SMC data in lGrimm et all 
(2003) wa s based on the star f ormation rate estimates ob¬ 
tained bv lFilioovic et all dl99SF) from the supernovae rate. 

5.3 The XLF shape 

At the first glance the shape of the HMXB XLF obtained 
by XMM-Newton appears to be more complex than a single 
slope N(> L) oc L~ 0 6 power law in the entire ~ 10 34 — 10 38 
erg/sec luminosity range. 

The maximum likelihood fit to the bright end of the 
distribution (Lx > 2 ■ 10 35 erg/s) with a power law results in 
the best fit value of the cumulative slope a = 0.74Aq jg which 
is consistent with the slope of 0.6. This is also confirmed by 
the Kolmogorov-Smirnov (K-S) test which gives probability 
of ~ 17%. Thus, the bright end of the XLF is consistent 
with the universal HMXB XLF. This conclusion is further 
confirmed by the ASCA XLF containing by a factor of ~ 3 
more sources. 

At low luminosities the uncertainties of HMXB identi¬ 
fication procedure become of importance and we consider 
the upper and lower limits of the XLF separately. The up¬ 
per limit on the XLF (the upper histogram in Fig. EJ is 
consistent with a single slope power law in the luminos- 
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ity range Lx > 10 34 erg/s with the K-S test probability 
of ~ 75% and the best fit value of the slope a = 0.5 ± 0.08. 
For the lower histogram in Fig. 0 (luminosity distribution 
of reliable HMXB candidates) we obtain best fit slope of 
a = 0.37 ± 0.08 in the full luminosity range. The observed 
distribution, is consistent with a single slope power law - the 
K-S test probability is « 33%. A power law approximation 
in the 10 34 — 2 ■ 10 35 erg/sec luminosity range gives the best 
fit slope of a O.ISIq’ 33 . 

Thus, the behavior of the XLF in the faint end is some¬ 
what uncertain. However, based on the statistical arguments 
and the distribution of the optical magnitudes (see sec. 14.31 
and 15.41 . the majority of sources of uncertain nature are 
expected to be the result of chance coincidence. Therefore 
the lower histogram in Fig. 0 should better represent the 
true HMXB XLF and the flattening of the XLF in the faint 
end can be r eal. It can result from the “prop eller” effect, as 
suggested by IShtvkovskiv fc Gilfanov! • 2 ll(l~Ji . 

5.4 Optical properties of HMXB candidates 

The distribution of V-band magnitudes of the HMXB can¬ 
didates is shown in Fig. |fi] In the case of the likely HMXB 
candidates it shows bimodal structure with two peaks at 
my ~ 14.6 and my « 15.5. To estimate the statistical 
significance of the observed bimodali ty, we perfo rm a test 
based on the bootstrap method ISilvermanlllQSlln . The ob¬ 
tained value of significance is rather low and does not allow 
us to make a definite conclusion (with probability that the 
observed distribution is multimodal of 97%, i.e. significance 
slightly above 2 <r). 

The origin of this bimodality, if real, is not clear. It 
could correspond to Be/X-ray systems with main sequence 
and giant Be secondaries. However, this explanation is not 
satisfactory as it would require a rather narrow distributions 
of the optical companions over spectral classes. 

The distribution of the my magnitudes of the uncertain 
HMXB candidates differs from those of likely candidates. 
This is consistent with expectations, that significant part of 
them should be a chance coincidences (see sec. 14.31 . 


6 LOG(N)—LOG(S) DISTRIBUTION OF CXB 

SOURCES IN THE DIRECTION OF SMC 

The log(N)-log(S) distributions of CXB sources, obtained 
after removal of likely HMXB candidates only and all HMXB 
candidates are shown in FigQ We fit the resulting log(N)- 
log(S) distribution in the flux range Fx > 2 • 10 -14 erg 
cnU 2 s ~ 4 with a power law model N(> S) = k(S/So)~ a , 
where So = 2 ■ 1CU 14 erg/s/cm 2 . Best fit value of slope for 
the case with only likely HMXB candidates excluded is a = 
1.48 ± 0.12. For the case with all the HMXB candidates 
excluded we obtain a = 1.55 ± 0.13. 

We also compare the resulting log(N)-log(S) distribu¬ 
tions in the flux range Fx > 2 • IIP 1 4 erg cm - 2 s ~ 4 with 
CXB number counts obtained bv iMoretti et all (l200Bh . The 
Kohnogorov-Smirnov test accepts their model giving prob¬ 
abilities ~ 96% for the case with likely HMXB candidates 
excluded and ~ 76% for the case with all the HMXB candi¬ 
dates excluded. The nu mber of sourc e s at flux 2 • lO -14 erg 
cnP 2 s ~ 4 predicted bv IMoretti et all ll2003lh N= I 78 I 4 ®, 
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S x , 2—10 keV flux, erg/s/cm 2 


Figure 7. Cumulative logN—logS distribution of CXB sources, 
obtained after removal of all HMXB candidates (lower histogram) 
and only likely HMXB candidates (upper histogram). The_solid 
li ne sh ows the CXB log(N)-log(S) distribution from lUoreit.i et. all 
foon-fl . The brightest source has flux Sx = 2.85-1CP 12 erg/s/cm 2 . 

also consistent with the observed numbers, N= 167±13 and 
N= 186 ± 14 for the two above cases respectively. 

Thus, the CXB log(N)-log(S) in the direction of the 
SMC is consistent both in the shape and the normalization 
with results of other CXB surveys. 


7 SUMMARY 

Based on the archival data of XMM-Newton observations, 
we studied the population of point-like sources in the field 
of the SMC. The total area of the survey is w 1.5 sq.degr. 
with the limiting sensitivity of ~ 10~ 14 erg/s/cm 2 (Fig|5J, 
corresponding to the luminosity of ~ 2.3 • 10 33 erg/s at the 
SMC distance. 

(i) Out of the 196 compact sources, detected in the 2-8 
keV energy band, ~ 3/4 are CXB sources, observed through 
the SMC (section 13. II Fig|3J. 

(ii) Based on the stellar mass and the star formation rate 
of the SMC we demonstrate, that the majority of the in¬ 
trinsic SMC sources, detected in the 2-8 keV band are high 
mass X-ray binaries (section . 

(iii) The proximity of the SMC and adequate angular res¬ 
olution of XMM-Newton make it possible to reliably fil¬ 
ter out the sources, whose properties are inconsistent with 
HMXB nature. Based on the optical and infrared magni¬ 
tudes and colors of the optical counterparts of the X-ray 
sources we identify 32 likely HMXB candidates (29 of which 
were previously known) and 18 sources of uncertain nature 
(section Table |5J. The remaining ~ 146 sources, with 
a few exceptions, are background objects, constituting the 
resolved part of CXB. Their flux distribution is consistent 
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both in shape and normalisation with other determinations 
of the CXB log(N)-log(S) (section [IJ Fig0. 

(iv) With these results we constrain lower and upper 
bounds of the luminosity distribution of HMXBs in the ob¬ 
served part of the SMC. We compare these with the extrap¬ 
olation towards low luminosities of the universal luminosity 
function of HMXBs, derived bv lOrimm et alJ ll2003fl (Fig.|SJ. 

The observed number of HMXBs is consistent with the 
prediction based on SFR estimates derived from supernovae 
frequency and analysis of the stellar color-magnitude dia¬ 
grams. If, on the contrary, the true value of the SFR is bet¬ 
ter represented by FIR, H a and UV based estimators, then 
the abundance of HMXBs in the SMC may significantly (by 
a factor of as much as ~10) exceed the value derived for the 
Milky Way and other nearby galaxies (sections 13.31 IHt. 

The shape of the observed distribution at the bright end 
is consistent with the N(> L) oc L -0 ' 6 power law of the 
universal HMXB XLF. At the faint end, Lx < 2 ■ 10 35 erg/s, 
the upper limit of the luminosity function is consistent with 
the L~ 0 ' 6 power law, while the lower limit is significantly 
flatter. 
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APPENDIX A: CROSS-CORELLATION OF 
TWO CATALOGUES 

Below we consider the problem of cross-corellation of two 
catalogs (e.g. of the list of detected X-ray sources and an 
optical catalogue). More specifically we concentrated on the 
case when only a fraction of X-ray sources has true optical 
counterparts in the optical catalogue and when the probabil¬ 
ity of chance coincidences is not negligible. In this case, the 
matches of X-ray sources with optical sources will include 
real optical counterparts and chance coincidences. We con¬ 
sider the behavior of the growth curves - the dependence of 
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the number of matches on the search radius. We show that 
using different dependences of the number of true matches 
and chance coincidences on the search radius, it is possible 
to estimate the number of true optical counterparts. 

A1 Case of constant density of optical stars 

The number of X-ray sources with optical matches inside ra¬ 
dius r depends on the number of X-ray sources which have 
true optical counterpart inside r and number of chance co¬ 
incidences of X-ray sources with field stars: 

Nxmatch = N( 1 - e~ pnr2 ) + Me - '”"' 2 [ <t>(r')dr', (Al) 

Jo 

where N - total number of X-ray sources, M - number of 
X-ray sources with true optical counterparts, p - density 
of optical field stars (i.e. excluding the true optical coun¬ 
terparts) and (j>{r) is the distribution of distances between 
X-ray sources and their optical counterparts. The first term 
in the above equation is a number of chance coincidences, 
i.e. number of X-ray sources which have at least one field 
star inside its search radius. The second term corresponds 
to number of true matches falling inside search radius r. The 
exponent in this term excludes cases when both true opti¬ 
cal counterpart and chance coincidence with field stars are 
present, as such matches were already taken into account 
in the first term. Assuming that all X-ray sources have the 
same positional error a (hence cf>(r) = r /a 2 • e~ r / 2<T ), we 
can rewrite this equation in the following way: 

N Xmatch = JV(1 - e~ pnr2 ) + Me-^(1 - e^'^ 2 ). (A2) 

Fitting the curve Nxmatch (r) with this formula one can 
determine the number of X-ray sources with true optical 
matches, surface density of optical stars p and error on the 
source positions. 

A2 Non-constant surface density of optical stars 

In a real situation, the density of optical stars p varies from 
field to field and even from source to source. This can lead 
to significant deviations of Nxmatch (r) curve behaviour from 
the eq. IA2I There is no obvious way to overcome this prob¬ 
lem completely. However, the density variations could be 
partially taken into account by replacing the mean density 
of optical objects in the first term of ea, IA2l bv the local den¬ 
sities of optical objects around each individual X-ray source. 
The first term in this case, of course, should be replaced by 
the corresponding sum over all X-ray sources. The density 
p in the second term, however, cannot be straightforwardly 
calculated because it is related to X-ray sources with true 
optical counterparts, which we do not know a priori. Thus, 
modifying the first term in equation I A2I we obtain: 

N 

N Xmatch = ^(l-e- w " 2 ) + Me-^ 2 (l - e — 2/2 - 2 ) . (A 3) 

i= 1 

The local densities pi of optical objects should be calculated 
individually for each X-ray source. The areas over which 
these densities are calculated must be large enough to min¬ 
imize the relative contribution of the true optical counter¬ 
parts and Poissonian errors. Density p corresponds now to 



Figure Al. The number of X-ray sources with an optical match 
as a function of search radius and its approximation with eq. ED 
The individual components of this equation corresponding to spu¬ 
rious matches with field stars and true optical counterparts (which 
don’t have simultaneous spurious matches) are shown with dashed 
and dot-dashed lines respectively. 


the mean density of optical objects near X-ray sources which 
have true optical counterparts. 

Fit parameters are the number of X-ray sources with 

true optical matches M, the value of positional error a and 

the mean density of optical objects p. Fitting the Nxmntek(r) 

curve using real list of XMM sources and optical sources 

from the MCPS catalogue we obtain the best fit value for 

the total number of sources with true optical counterpart 

in our sample is ~ 31, surprisingly close to the number of 

our likely HMXB candidates (see Fig. ED . Best fit value of 

positional error, a ss 1.2", also agrees with typical values for 

X-ray sources from our list ranging from « 0" to « 2". We 

note also, that from the eq. [H it follows that subtracting the 

first sum from Nxmatch curve and multiplying the remaining 
2 

part by e pnr , we can obtain the dependence of the number 
of true counterparts among optical matches on the search 
radius. The resulting curve is shown on Fig. EH It is clear 
from this figure, that almost all true optical counterparts of 
X-ray sources will be detected with the search radius of 4" 
used in our HMXB identification procedure. 


A3 Account for spread in distribution of 
positional uncertainties 

The deviation of positional errors distribution from the delta 
function together with the astrometric uncertainty can also 
modify Nxmnfck(r) curve. With the spread in positional er¬ 
rors accounted for, the equation 1X51 becomes: 

N Xmateh = £)(1 - e~ Pinr2 ) + Me""" 2 • (A4) 
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Figure A2. Number of true matches as a function of radius 
estimated using eq. EH (see E2J- The vertical dashed line shows 
the search radius of 4 arcsec. It is clear from this figure, that 
almost all true optical counterparts of X-ray sources fall within 
the chosen search radius. 



<r)g(cr)dcrdr ', 


_ 2 /<-> 2 

where a) oc re~ r ' a is the (Gaussian) distribution of 
distances between the X-ray and optical sources for given a 
and g(c r) is a distribution of positional uncertainties of X- 
ray sources. To show that such a spread will not significantly 
change our results, we assumed that errors distribution fol¬ 
lows the power law g(a) oc cr“ between cr = 0.2” and 2.3”. 
In this case the free fit parameters are the number of X-ray 
sources with true counterparts M, optical sources density p 
and slope a. 

Fitting our curve with this formula we obtained the best 
fit value M=34, which also agrees well with the number of 
our likely HMXB candidates. 









